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Abstract 

The unique properties of central potential of the form —f3e~ r r 1 were studied 
using the recently developed critical parameter technique. The particular cases 
of 7 = and 7 = — 1 yield, respectively, the exponential and Yukawa potentials 
widely used in the atomic, molecular and nuclear physics. We found different 
behavior of the energy levels of this potential for three different ranges of the 
value of 7. For 7 > it was found that the energy of bound states with the same 
principal quantum number N decreases with increasing angular momentum t. 
The Gaussian and Woods-Saxon potentials also show this behavior. On the 
contrary, for — 2 < 7 < — 1 increasing i gives a higher energy, resembling the 
Hulthen potential. However, a potential with — 1 < 7 < possesses mixed 
properties, which give rise to several interesting results. For one, the order of 
energy levels with different quantum numbers is not preserved when varying the 
parameter /3. This leads to a quantum degeneracy of the states, and in fact, 
for a given value of 7 we can find the values ftthr for which two energy levels 
with different quantum numbers coincide. Another interesting phenomena is 
the possibility, for some values of 7 in this range, for two new energy levels 
with different quantum numbers to appear simultaneously when f3 reaches their 
common critical value. 

Keywords: central potentials, critical parameters, energy levels, quantum 
numbers, angular momentum 



1. Introduction 

The interest in the ordering of energy levels is as old as quantum mechanics. 
It can be reflected in the electron configurations of the elements in the periodic 
table, and in the nuclear shell model. In the 1980s there was a renewed interest in 
this subject following the discoveries of mesons composed of the heavier quarks 
c and b. From the disposition of the energy levels of these mesons, the nature of 
the quark-quark (or, rather quark-antiquark) potential can be deduced, if one 
can relate the properties of the potential with the ordering of the energy levels 



* Corresponding author 
Email address: livertsOphys.huji.ac.il (Evgeny Z. Liverts) 



Preprint submitted to Elsevier 



March 4, 2013 



pertaining to it. This development led the authors of Q (following and 0]) 
to prove that for non-relativistic two-body systems the order of energy levels 
with the same principal quantum number N is controlled by the sign of the 
Laplacian of the spherically symmetric potential V(r), 

. , N 2dV d 2 V 

where N — n + £ + 1, n is the number of nodes in the radial wave function 
and £ is the angular momentum quantum number. Specifically, they show that 
under the following two conditions the ordering of energy levels with the same 
principal quantum number N can be exactly determined^ 

1. If 

AV(r) > for all r > (2) 

then 



E n f > E n —H+l 



(3) 



2. If 



AV(r) < for r < r , and — < for r > r > 0, (4) 

dr 



then 



E n ,i < En-ij+i (5) 

These results were also generalized for the relativistic case in Ref. Q, but in 
this letter we limit ourselves to the non-relativistic realm. 

The question arises from these results - what can one say about the ordering 
of energy levels relating to potentials that do not fulfill either of the conditions 
© or (HI above? In this letter we try to answer this question for the exponential- 
power (EP) potential, i.e., of the form 

V lh (r) = -pe- r Ti G8>0, 7 >-2). (6) 

This is done using the critical parameter technique (CPT) recently introduced 
in where it was also used to study the energy levels of a few short-range 
central potentials widely used in atomic, molecular and nuclear physics. Two 
of these, namely the exponential and Yukawa potentials are special cases of Va 7 
with 7 equals and —1, respectively. 



1 Notice that for the 2 nd case only the weaker condition is given here, instead of the more 
restrictive condition, AV(r) < for all r > 0, which is the opposite of (2~l . 
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2. Methods 



In [f| the CPT was introduced as a methodology for obtaining information 
on the bound solutions of the non-relativistic Schrodinger equation with central 
potentials of the form 

V (r) = -V f (r/r s ) , (V > 0) (7) 

where r s is the so called screening parameter, and Vq represents the coupling 
constant. Taking into account the scaling properties, the solution of the relevant 
equation depends effectively on a single parameter /3 = 2mgr 2+1 / ti 2 , where m is 
the reduced mass of the considered system. Parameter g, which is proportional 
to Vo, will be defined in ©. The CPT consists primarily of finding at what 
values of /? = Pn,i, a s it is being varied, an eigenstate of the Hamiltonian 
with given {n,C\ becomes a transition state, i.e., will have zero energy. Further 
increasing j3 beyond this "critical" value (CP) will make this eigenstate a bound 
state, causing a new bound energy level E n ^ to "appear" . 

This technique allows one to answer such questions as a) How many bound 
states with given I exist for a given potential? b) What is the maximum I a 
bound state can have, given the potential? etc. 

In [f| it is also shown that the CPT can determine the order in which the 
different energy levels "appear", as j3 is increased. This led the authors to 
conjecture that this order of the energy levels is preserved when /3 is further in- 
creased beyond the values at which the levels first appeared. If this conjecture 
is true, then one can very easily determine the ordering of energy levels of ever 
potential that can be solved using the CPT. For example, the authors of 
derive the conditions (67), regarding the ordering of the energy levels of the po- 
tentials studied there, by relying only on their knowledge regarding the relations 
between the CP's obtained for each of these potentials. Our results suggest that 
the validity of the conjecture depends on the behavior of the potential near the 
origin, as shown in several examples below. 

A few technical details are due here. All calculations were provided by 
the Mathematica 8 Q codes which realized methods described in [&]. The 
CP's depend on the asymptotic behavior (r —> oo) of x( r X the solution to the 
corresponding Schrodinger equation 



d 2 x(r) 
dr 2 

when 



-I3v(r) + - E nA (J3) 



E n , e 08) — — > 0. 



X(r) (8) 



However, the numerical integration of the above differential equation starts from 
some small value of r. Therefore, the accuracy of the asymptotic behavior is 
defined by our possibility to calculate accurately the radial wave function (and 
its derivative) at this initial value of r which should be chosen as far as possible 
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from the origin. Let the leading term of a series expansion for the potential 
under consideration has the form: 



v(r) ^ A- (5>0, a<2). (9) 

A series expansion for the reduced radial wave function in the case of a potential 
with integer or half- integer a was presented in In order to provide calcula- 
tions of the CP's, e.g., for the EP potential with any real 7 = —a > —2, one 
needs a series expansion for the corresponding radial wave function. Therefore, 
we present such an expansion: 



C (r) = Cr^ 1 + iP)r k{2 - a)+l . (10) 



.1$ ( flV k (2"«)+i 
\ k=l i=0 / 

Here C is an arbitrary constant, and the coefficients A^(/3) represent polyno- 
mials in /3. The duplicated powers of r certainly have to be dropped. Sub- 
stituting Eq. (|TU1) into Eq. ©, and then equating sequentially (starting from 
the lowest power) the expansion coefficients of the same powers of r for the 
left-hand and right-hand sides, one can calculate any finite number of the co- 
efficients A$(0). This expansion can be used to solve Eq.® with or without 
taking the limit E n< £ —> in order to find either the CP's or the energies £"„^(/3), 
respectively. 

It is interesting to note that even though the leading term of this expansion, 
r f+1 , is well-known, we could not find the general form (fTOj) in the scientific 
literature. 

3. Cases 

It is easy to check that the Hulthen and the Yukawa potentials fulfill the 
condition (0} , and therefore their energy levels should follow ([3]) . This coincides 
with the results presented in [B| (see the first condition in (67) there). 

On the contrary, the energy levels of the exponential, the Gaussian and 
the Woods-Saxo potentials should satisfy the inequality Q , according to the 
second of conditions (67) in [5]. Interestingly, however, the Laplacians and the 
first derivatives of these potentials obey the conditions 

AV(r) % for r ^ r , 
dV(r)/dr > for all r > 1 ' 

which differ slightly from ([2]), and in fact do not coincide with any of the con- 
ditions mentioned in [l[. Our results for these potentials satisfy inequality ([3]), 
in accord with the conjecture. 



2 However, most of the coefficients are monomials. 

3 Two versions of the Woods-Saxon potential with different parameters were studied in [j| 
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Note that the Yukawa and exponential potentials 0, which are the special 
cases with 7 = —1, of the EP potential ([6]), demonstrate opposite ordering for 
energy levels with equal principal quantum numbers. The Laplacian of potential 
([B]) has a form0: 



AV(r) = V(r) 



( 7 +l)(7-2r) 



(12) 



It is easily seen that for 7 < — 1 the right-hand side of Eq. (IT^I) is negative for all 
r > 0. Hence, the corresponding energy levels must satisfy inequality ([5]). On 
the other hand, it follows from the results of [H[ that for 7 = 00 the inequality ([3]) 
holds. This is connected to the fact that the leading term in expansion of the 
exponential potential near the origin is proportional to zero or positive power of 
r (such is also the case for the Gaussian and Woods-Saxon potentials, mentioned 
above). This behavior suggests that the potential © with — 1 < 7 < may 
possess some mixed properties. 

In order to check this we have computed the CP's for the EP potential ([6]) 
with different 7 in the range (—1, 0). 



4. A counterexample 

The CP's for EP-potential with 7 = —1/2 are presented in Table Q] One 
can observe that the critical parameters (3 n .i and fi n -x,i+i (that relate to the 
same principal quantum number TV) are on the diagonals connecting the {n, £}- 
positions assigned as {ro,0} on the lower left and {l,n — 1} on the upper right. 
The smallest parameter /3 m i n of each diagonal was framed. 

The critical parameters determine the order in which the energy levels ap- 
pear. According to the conjecture formulated in section 2, it is "naturally" 
assumed that this order is preserved when the parameter f3 of the EP poten- 
tial ([6]) is increased. 

If this assumption is true, the energy levels E n> t (/3) corresponding to (3 > 
@ n ^i must satisfy the inequality <j5j) for {n, £} belonging to CP's which are framed 
or are above the framed ones. On the contrary, the energy levels for the CP's 
which are below the framed ones, must obey the inequality (|3]). 

The computational results presented in Table [2] show that this assumption 
does not hold here. It follows from Table [Tj that for the EP potential ([6]) with 
parameters 7 = —1/2 and j3 > /3i t g = 194.393 the Schrodinger equation ([8]) has 



4 It is worth noting that the potential of the form (0 was also studied by means of the so 
called auxiliary field method in [f| . There the critical parameters (or "critical heights" in the 
language of |(|) were obtained for the exponential and the Yukawa potentials with £,n < 3, 
and within the limits of their inferior precision, coincide with the results from |5|. 

5 Note that the Laplacian of the Yukawa potential (7 = — 1) reduces to the Yukawa potential 
itself. 

6 It is easy to guess that for 7 > 0, as well. 

7 It is seen that the smallest parameters occupy the positions {n, 1} with n = 21 + p, where 
V = 1,2,3. 
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bound state solutions only for £ < 8 and n < (9 — £). The corresponding binding 
energies E n ^i presented in Table [2] were computed by two methods for fidelity. 
First, the results were obtained bysolving Eq. ([5]) directly. Then, application 
of the technique presented in Ref. [8J enabled us to verify the accuracy of these 
results. Table [5] demonstrates that only the four energy levels £" n /(200) with 
n + £ — 9 and 5 < £ < 8 satisfy the inequality (JS|), whereas levels with 
< I < 5 (and the same N = 10) fulfill the opposite condition ([3]). This 
contradicts the assumption, and disproves the conjecture. 

Where does it fail? For the conjecture to hold, the energy levels must re- 
tain the order in which they appear, as the universal parameter /? is increased. 
Obviously, all existing energy levels are lowered when (3 increases, until a new 
level appears. However, for the EP potential (O with —1 < 7 < the "speed of 
lowering" the energy levels E n i and E n >^> with S = n — n' = £' — £ — ±1 (and 
the same TV), is opposite to the order of their appearance in this process. That 
is, a level E n ^ that has just appeared, will be going down in energy "faster" 
than the originally lower level (with the same N and 6 = ±1) that had appeared 

(n £) 

before it. Therefore, at some value of j3 = pl^ r these levels will reverse their 
order. 



5. Degeneracy and simultaneous appearance 

The fact that the "speed of lowering energy" of the levels that had just 
appeared is higher than that of levels with the same ./V and 6 = ±1 but that 
had appeared before hand, has several interesting consequences. 

First of all, as suggested above, this behavior will cause a degeneracy of 
two energy levels E n j^ — E n ^u + i, when f3 reaches the value for which these 
two levels reverse their order. Indeed, for the example above (with 7 = —1/2 

and j3 > /3i,s), at some threshold /3 = tne energy E^s (Pthr^ ke- 

comes equal to _E 3j6 (p^j • The further increase of /3 to P^'P > ftfyj^ 

leads to E 2 , 7 (X ( '; 7) ) = £3,6 (til'?) < E 4 . 5 (^ t ( '; 7) ). And at last, for some 

& = & > & one obtains E hS (/3^ 8) ) = #2,7 (til?) < #3,6 (&)■ 
(1 8) 

Finally, for f3\^ r < j3 < /3i,g all possible energy levels with N = 1Q satisfy 
the condition The situation is of course similar for smaller j3. The corre- 
sponding Pthr are presented in Tabled where £ ma x is the largest orbital number 
which admits a bound state for the given parameter /3 — fithr- It is seen that 
for N= £ max + n = 8, 9 there are 3 threshold potentials (n = 1, 2, 3) for each N. 
For N= 5, 6, 7 there are only 2 threshold potentials (n = 1, 2) for each N. And 
at last, for N= 2,3,4 one obtains the threshold potentials producing the pairs 
of states with _ 

E 3p (til iy ) = EzM'A Ea4&) = EaMh?) andS 5/ (/3 t t 3) ) = E, d {&) 
with the jsame principal quantum numbers. Here we used the widespread nota- 
tion for E N ,i = E n j. 
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The second consequence of this behavior is the possibility for simultaneous 
appearance of two levels, i.e., there are certain values of — 1 < 7 < for which 
two new levels E n j = -EVs-i, £+i(~ —0) "appear" simultaneously at the same 
value fi n .g — fi n -i,t+i- Our computations show that, e.g., 

#2,o = - 7-9797 for 7 = -0.238825, (13) 
£3,0 = #2,1 ^ 19.1Q36 for 7 = -0.42046, (14) 

and so on. 

How does this behavior depend on the parameters of the potential ©? Recall 
that the framed values in Table Q] represent the smallest for all levels with the 
same N and 6 = ±1. For 7 = — 1/2 they are situated along the "diagonal" curve 
of the table. As 7 approaches 0, the curve of framed (smallest) (3 approaches the 
lower left corner. At the limit 7 = all the potentials will satisfy ([3]), regardless 
of j3. On the contrary, as 7 approaches —1, the curve of smallest (3 approaches 
the upper right corner, and at the limit 7 = —1 all the potentials will satisfy §5§. 

Also note that beyond these limits (7 < — 1, or 7 > ) the "speed of lowering 
energy" of the levels with the same N and different £ now matches the order of 
their appearance. This is related to the behavior of the potential near the origin. 
As was mentioned above, this is the reason why the exponential, Gaussian and 
Woods-Saxon potentials satisfy ([3]), even though they only fulfill the conditions 
(fTT]> rather than ©. 

6. Conclusions 

We have used the CPT developed in [5j , together with a new series expansion 
for the radial wave function (|10[) , to study the ordering of energy levels with the 
same principal quantum number N for a central potential of the form ([5]), as 
well as for several other short-range central potentials widely used in different 
fields of physics. We found that our results are not only consistent with the 
conditions formulated in [l[, but they even suggest a possible generalization of 
the condition ([2]) to (fTTj) . We found that the ordering of levels of the potential 
([6]) can be characterized according to three ranges of the values of the parameter 
7 

E n ,e > E n _i /+ i for 7 > 

E n< i = E n _ 1>t+1 for - 1 < 7 < 
E ni£ < £„_ M+ i for - 2 < 7 < -1 

Although it was not shown explicitly, these different properties are related 
to the behavior of the potential near the origin, and perhaps this rule can also 
be generalized to other potentials. 

In the range — 1 < 7 < the above properties give rise to several interesting 
consequences for the EP potential ([6]), such as the degeneracy of two (or more) 
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energy levels E n< i (/3) = E n -x,l+\ (P) at some calculable values of /? = /Sthn or 
the simultaneous appearance of two new energy levels (as transition states, cf. 
[5]) at the same value of the critical parameter j3 = fi n ,i — ftn-i.e+i- 
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Table 1: Critical parameters /3„ p of the EP-potential (O with 7 = —1/2, V(r) = 
—/3exp(—r)/y/r. 



n\£ 





1 


2 


3 


4 


5 


6 


7 


8 


9 


1 


1.72515 


8.77711 


20.6738 


37.4369 


59.0736 


85.5867 


116.977 


153 


246 


194.393 


240.419 


2 


7.95958 


19.2665 


35.3988 


56.3915 


82.2555 


112.995 


148.611 


189 


105 


234.478 


284.729 


3 


18.8289 


34.2285 


54.5185 


79.6946 


109.755 


144.699 


184.525 


229 


232 


278.82 


333.289 


4 


34.341 


53.7313 


78.1078 


107.418 


141.639 


180.761 


224.775 


273 


678 


327.468 


386.142 


5 


54.498 


77.8093 


106.213 


139.61 


177.956 


221.226 


269.405 


322 


485 


380.46 


443.325 


6 


79.3008 


106.482 


138.864 


176.307 


218.742 


266.131 


318.451 


375 


687 


437.829 


504.872 


7 


108.75 


139.762 


176.081 


217.532 


264.024 


315.505 


371.942 


433 


313 


499.606 


570.809 


8 


142.845 


177.657 


217.879 


263.306 


313.825 


369.371 


429.901 


495 


388 


565.812 


641.161 


9 


181.588 


220.174 


264.269 


313.642 


368.162 


427.749 


492.35 


561 


931 


636.47 


715.947 


10 


224.977 


267.316 


315.259 


368.553 


427.049 


490.653 


559.304 


632 


961 


711.596 


795.187 


11 


273.013 


319.087 


370.857 


428.048 


490.497 


558.097 


630.778 


708 


493 


791.207 


878.895 


12 


325.696 


375.489 


431.066 


492.134 


558.515 


630.091 


706.784 


788 


538 


875.315 


967.086 


13 


383.027 


436.525 


495.891 


560.818 


631.112 


706.646 


787.332 


873.11 


963.934 


1059.77 


14 


445.004 


502.195 


565.336 


634.104 


708.294 


787.769 


872.432 


962.217 


1057.07 


1156.97 


15 


511.629 


572.501 


639.403 


711.998 


790.068 


873.467 


962.091 


1055.87 


1154.74 


1258.67 


16 


582.901 


647.445 


718.096 


794.502 


876.437 


963.745 


1056.32 


1154.07 


1256.95 


1364.91 


17 


658.82 


727.027 


801.414 


881.62 


967.407 


1058.61 


1155.11 


1256.83 


1363.7 


1475.68 


18 


739.386 


811.249 


889.362 


973.355 


1062.98 


1158.07 


1258.49 


1364.16 


1475.01 


1590.98 


19 


824.6 


900.11 


981.94 


1069.71 


1163.16 


1262.12 


1366.45 


1476.06 


1590.87 


1710.84 


20 


914.46 


993.611 


1079.15 


1170.68 


1267.95 


1370.77 


1478.99 


1592.53 


1711.3 


1835.24 



Table 2: Binding energies — E nt i for the EP-potential V(r) = — 200 exp(— r)/y/r. 



n\£ 





1 


2 


3 


4 


5 


6 


7 8 


1 


453.076 


253.517 


161.668 


106.560 


69.5605 


43.3235 


24.2864 


10.5122 


0.943082 


2 


217.692 


142.656 


95.1156 


62.4832 


39.1216 


22.1330 


9.89629 


1.52600 




3 


123.435 


83.5646 


55.2885 


34.7710 


19.7922 


9.04198 


1.79035 




4 


71.9016 


48.0047 


30.3073 


17.3037 


7.99882 


1.81080 






5 


40.6508 


25.7661 


14.7085 


6.81296 


1.64754 








6 


21.1813 


12.0487 


5.53120 


1.35410 










7 


9.36962 


4.20250 


0.982547 












8 


2.88204 


0.58778 














9 


0.234378 
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Table 3: Threshold parameters f) thr of EP-potential with 7 = — 1/2 for E( n z ( max ) {Pthr) = 

E {n+l,t max -l)(Pthr)- 



n\{£max + n) 


2 


3 


4 


5 


6 


7 


8 


9 


1 


9.002821 


22.22577 


41.3813 


66.4516 


97.4325 


134.3227 


177.1216 


225.82855 


2 








59.9317 


89.30232 


124.6208 


165.8688 


213.0374 


3 














155.5733 


201.1049 
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